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We perform the first studies of various inter-quark potentials in SU(3) C lattice QCD. 
From the accurate lattice calculation for more than 300 different patterns of three-quark 
(3Q) systems, we find that the static 3Q potential V3Q is well described by Y-Ansatz, i.e., 
the Coulomb plus Y-type linear potential. Quark confinement mechanism in baryons is 
also investigated in maximally-Abelian projected QCD. We next study the multi-quark 
potentials V„q (n=4,5) in SU(3) C lattice QCD, and find that they are well described 
by the one-gluon-exchange Coulomb plus multi-Y type linear potential, which supports 
the flux-tube picture even for the multi-quarks. Finally, we study the heavy-heavy-light 
quark (QQq) potential both in lattice QCD and in a lattice-QCD-based quark model. 
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1. Introduction 

In 1966, Yoichiro NambrP first proposed the SU(3) C gauge theory, i.e., quantum 
chromodynamics (QCD), as a candidate for the fundamental theory of the strong 
interaction, just after the introduction of the "new" quantum number, "color" !^ 
In 1973, the asymptotic freedom of QCD was theoretically provenj^ and, through 
the applicability check of perturbative QCD to high-energy hadron reactions, QCD 
has been established as the fundamental theory of the strong interaction. However, 
in spite of its simple form, QCD creates thousands of hadrons a nd le ads to various 
interesting nonperturbative phenomena such as color confinementPEl and dynamical 
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chiral-symmetry breaking.^ Even at present, it is very difficult to deal with QCD 
analytically due to its strong-coupling nature in the infrared region. Instead, lattice 
QCD has been applied as a direct numerical analysis for nonperturbativc QCD. 

In 1979, the first application of lattice QCD Monte Carlo simulations^ was done 
by M. Creutz for the inter-quark potential between a quark and an antiquark using 
the Wilson loop. Since then, the study of inter-quark forces has been one of the 
important issues in lattice QCD.^ Actually, in hadron physics, the inter-quark force 
can be regarded as an elementary quantity to connect the "quark world" to the 
"hadron world" , and plays an important role to hadron properties. 

In this paper, we perform the first detailed studies of the various inter-q uark 
forces in the three-quark and the multi-quark systems with SU(3) C lattice QCD.^UHl 

2. The Three-Quark Potential in SU(3) Lattice QCD 

In general, the three-body force is regarded as a residual interaction in most fields 
in physics. In QCD, however, the three-body force among three quarks is a "pri- 
mary" force reflecting the SU(3) gauge symmetry. In fact, the three-quark (3Q) 
potential is directly responsible for the structure and properties of baryons, similar 
to the relevant role of the QQ potential for meson properties. Furthermore, the 3Q 
potential is the key quantity to clarify the quark confinement in baryons. However, 
in contrast to the QQ potential,^ there were almost no lattice QCD studies for the 
3Q potential before our study in 1999,^ in spite of its importance in hadron physics. 

As for the functional form of the inter-quark potential, we note two theoretical 
arguments in short and long distance limits. 

1. At short distances, perturbative QCD is applicable, and the inter-quark potential 
is expressed as a sum of two-body one-gluon-exchange (OGE) Coulomb potentials. 

2. At long distances, the strong-coupling expansion of QCD is plausible, and it leads 
to the flux-tube picture. 

Then, we theoretically conjecture the functional form of the inter-quark potential as 
the sum of OGE Coulomb potentials and the linear potential based on the flux-tube 
picture. Of course, it is highly nontrivial that these simple arguments on UV and 
IR limits of QCD hold for the intermedia te re gion. Nevertheless, the QQ potential 
Vqq(r) is well described with this form as^^l 

^qqW = -^ + W + c qq- (!) 

For the 3Q system, there appears a junction which connects the three flux-tubes 
from the three quarks, and Y-type flux-tube formation is expected.^ Then, the 3Q 
potential is expected to be the Coulomb plus Y-type linear potential, i.e., Y-Ansatz, 

V3Q = ~^3Q ^ 1 Z 1 + ^Q^mm + C 3 Q, (2) 

i<j < ■>' 

where L m i n is the minimal total length of the Y-shaped flux-tube. 
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2.1. Manifestly Gauge-Invariant Lattice QCD Results 

For more than 300 different patterns of spatially-fixed 3Q systems, we calc ulate 
the 3Q potential V 3Q from the 3Q Wilson loop W SQ using SU(3) lattice QCE0E3 
with the standard plaquette action at the quenched level on various lattices, i.e., 
09=5.7, 12 3 x 24), (/3=5.8, 16 3 x 32), (/3=6.0, 16 3 x 32) and (J3 = 6.2, 24 4 ). For the 
accurate measurement, we construct the ground-state-dominant 3Q operator using 
the smearing method. Note also that the lattice QCD calculation is completely 
independent of any Ansatz for the potential form. 

To conclude, we find that the static 3Q potential V-jq is well described by the 
Coulomb plus Y-type linear potential (Y-Ansatz) within 1%-level deviation!^ We 
also find the universality of the string tension, CT3Q ~ ctqq, and the OGE result, 
A^q ~ Aqq/2. As an example, we show in Fig. 1(a) the 3Q confinement potential 
^3Q nf > i- e -' the 3Q potential subtracted by the Coulomb part, plotted against the 
Y-shaped flux-tube length L m i n . At each (3, clear linear correspondence is found 
between V 3 Q nf and L nlin , which indicates Y-Ansatz for the 3Q potential.^ 



3Q Confinement Potential [a" 1 ] 
(3Q Potential - Coulomb part) 




Lmin [a] 




Fig. 1. 



(a) (Left) The 3Q confinement potential V 3 C Q nf , 



i.e., the 3Q potential subtracted by the 



Coulomb part, plotted against the total flux-tube length L m % n of Y-Ansatz in the lattice unit. 10 (b) 
(Middle and right) Lattice results for Y-type flux-tube formation in the spatially-fixed 3Q system 
in MA projected QCD. The distance between the junction and each quark is about 0.5ml. 11 



2.2. Quark Confinement Mechanism in Baryons 

As another clear evidence for Y-Ansatz, Y-type flux-tube formation is actually 
observed in maximally- Abelian (MA) projected lattice QCD from the measurement 
of the action density in the spatially- fixed 3Q systempJJ as shown in Fig. 1(b). 

MA projected QCD includes electric and magnetic currents, which can be sepa- 
rated as the photon part and the monopole part with the Hodge decomposition [HEl 
We investigate these parts of the 3Q potential in MA projected QCD in quenched 
SU(3) lattice QCD with 16 4 and (3—6.0, and find the following result. 

1. The monopole part V 3 q° of the 3Q potential is almost single- valued function of 
L m i„, and is approximated as V^q° ~ 03QL m i n at long distances. 

2. The photon part V^q of the 3Q potential is almost single- valued function of 
Lcoui = (- + + -) — 1 with a, b and c being the three sides of the 3Q triangle. 
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Fig. 2. The lattice results for Hodge decomposition of the 3Q potential in MA projected QCD: 
(a) The monopole part V^q 1 plotted against L mln ; (b) the photon part V^q plotted against Lcoul- 



Thus, whereas the electric current induces Coulomb forces, the magnetic-monopole 
current induces the Y-type linear confinement potential for quarks in baryons. 

In this way, together with recent several analytical studies^U^ and other recent 
lattice QCD studies, Y-Ansatz for the static 3Q potential is almost settled. 



3. Inter-Quark Potentials in Multi-Quark Systems in Lattice QCD 

Next, we perform the first s tudy of the inter-quark interaction in multi-quark sys- 
tems in SU(3) lattice QCD.LIQ As for the potential at short distances, the per- 
turbative OGE potential would be appropriate, due to the asymptotic nature of 
QCD. For the long-range part, however, there appears the confinement potential as 
a typical nonperturbative property of QCD, and its form is highly nontrivial in the 
multi-quark system. In fact, to clarify the confinement force in multi-quark systems 
is one of the essential points for the construction of the QCD-based quark-model 
Hamiltonian. Then, we investigate the multi-quark potential in lattice QCD, with 
paying attention to the confinement force in multi-quark hadrons. 



3.1. One- Gluon- Exchange (OGE) Coulomb plus Multi-Y Ansatz 

We first consider the theoretical form of the multi-quark potential, since we will 
have to analyze the lattice QCD data by comparing them with some theoretical 
Ansatz. By generalizing the lattice QCD result of Y-Ansatz for the three-quark 
potential, we propose the OGE Coulomb plus multi-Y AnsatzjE3 

V nQ = 7^ A nQ 2_, 1 ' " I + ^"Q^min + C n Q (fl = 3, 4, 5, ...), (3) 

for the potential form of the multi-quark system. Here, the confinement potential 
is proportional to the minimal total length L m j n of the color flux tube linking the 
quarks, which is multi-Y shaped in most cases. 

In the following, we study the inter-quark interaction in multi-quark systems in 
lattice QCD, and compare the lattice data with the theoretical Ansatz ©. Note that 
the lattice QCD data themselves are meaningful as primary data on the multi-quark 
system directly based on QCD, and do not depend on any theoretical Ansatz. 
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3.2. Formalism of the Multi-Quark Wilson Loop 



We formulate the multi-quark Wilson loop to obtain the multi-quark potential in 
lattice QCD P^S Similar to the derivation of the QQ potential from the Wilson loop, 
the static multi-quark potential can be derived from the corresponding multi-quark 
Wilson loop. We construct the tetraquark Wilson loop W 4 q and the pentaquark 
Wilson loop W$q in a gauge invariant way as shown in Fig. 3, and define them as 



W 4Q = itr(Mi7?i 2 M 2 ii 2 ), 

W 5Q = l e abc e a ' b ' c ' M aa ' (R 3 R 12 R 4 ) bb ' (L 3 L 12 L 4 ) CC ' , 
where M, Mi, Lj, Rj (i=l,2, j=l,2,3,4), R 12 and L 12 are given by 
M, Mi , Rj , L 3 = Pexp{ig [ dx^A^x)} 6 SU(3) C 

JM, Mi, Rj, L, 



£>a'a ^abc^a'b'c' r>bb f r>cc' 

H 12 =2 1 2 ' 



j a 1 a ,afcc ,a' b' c' j bb' j cc' 



(4) 
(5) 

(6) 
(7) 



The multi-quark Wilson loop physically means that a gauge-invariant multi-quark 
state is generated at t = and annihilated at t — T with quarks being spatially 
fixed in R 3 for < t < T . The multi-quark potential is obtained from the vacuum 
expectation value of the multi-quark Wilson loop as 



^4Q = - lim ^ln(W 4Q ), 

1 — >OQ 1 



T iimIm(^ Q ) 



(8) 
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Fig. 3. (a) The tetraquark Wilson loop W4Q. (b) The pentaquark Wilson loop W5Q. The contours 
M, M{, Rj , Lj(i = 1,2, j = 3,4) are line-like and Rj,Lj(j = 1,2) are staple-like. 



3.3. Tetraquark Potential and Flip-Flop in SU(3) Lattice QCD 

For about 200 different patterns of QQ-QQ configurations, i.e., tetraquark (4Q) 
systems, we perform the detailed study of the tetraquark (4Q) potential Vaq in 
SU(3) lattice QCD with /?=6.0 and 16 3 x 32, and find the following results.^ 

1. When QQ and QQ are well separated, the 4Q potential V4Q is well described 
by the OGE Coulomb plus multi-Y Ansatz, which indicates the formation of the 
multi- Y-shaped flux-tube connecting four (anti)quarks as shown in Fig. 4(a). 

2. When the nearest quark-antiquark pair is spatially close, V4Q is well described by 
the sum of two QQ potentials, which indicates a "two-meson" state as Fig. 4(b). 
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Fig. 4. (a) A connected tetraquark (QQ-QQ) configuration and (b) a "two-meson" configuration, 
(c) A QQ-Q-QQ type pentaquark configuration. The lattice QCD results indicate the multi-Y- 
shaped flux-tube formation for the connected 4Q system and the QQ-Q-QQ system. 



As the examples, we show in Fig. 5 the lattice QCD results of the 4Q potential 
V4Q for symmetric planar 4Q configurations as shown in Figs. 4(a) and (b), where 
each 4Q system is labeled by d = Q1Q2/2 and h = Q1Q3. For large values of h 
compared with d, the lattice data obey the OGE Coulomb plus multi-Y Ansatz. 
For small h, the lattice data obey the "two-meson" Ansatz, where the 4Q potential 
is described by the sum of two QQ potentials, Vqq(ri3) + Vqq^a) = 2Vqq(/i). 

Thus, the tetraquark potential V4Q is found to take the smaller energy of the 
connected 4Q state or the two-meson stateP^ In other words, we observe a clear 
lattice QCD evidence of the "flip-flop", i.e., the flux-tube recombination between 
the connected 4Q state and the two-meson state. 
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Fig. 5. Lattice QCD results of the tetraquark potential V4Q for symmetric planar 4Q configu- 
rations in the lattice unit. The symbols denote the lattice QCD data. The solid curve denotes the 
OGE plus multi-Y Ansatz, and the dotted-dashed curve the two-meson Ansatz. 



3.4. Lattice QCD Result of the Pentaquark Potential 

We perform the first study of the pentaquark (5Q) potential V5Q in lattice QCD 
with /3=6.0 and 16 3 x 32 for 56 different patterns of QQ-Q-QQ type pentaquark 
configurations!^^ We find that the lattice QCD data of V5Q are well described by 
the OGE Coulomb plus multi-Y Ansatz, i.e., the sum of the OGE Coulomb term 
and the multi-Y-type linear term based on the flux-tube picture. 

We show in Fig. 6 the lattice QCD results of the 5Q potential Vsq for symmetric 
planar 5Q configurations as shown in Fig. 4(c), where each 5Q system is labeled 
by d = Q1Q2/2 and h = Q1Q3. In Fig. 6, we add the theoretical curves of the 
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OGE Coulomb plus multi-Y Ansatz, where the coefficients (^sQjOsq) are set to be 
(A 3 q,(T3q) obtained from the 3Q potentialPln Fig. 6, one finds a good agreement 
between the lattice QCD data of Vsq and the theoretical curves. 

In this way, the multi-quark potentials V n Q (n =3,4,5 ) are found to be well de- 
scribed by the OGE Coulomb plus multi-Y AnsatzP^Ell These lattice QCD results 
support the flux-tube picture even for the multi-quark systems. 




u T v 



Fig. 6. Lattice QCD results of the pentaquark potential V5Q for symmetric planar 5Q config- 
urations as shown in Fig. 4(c) in the lattice unit: (a) V 5 q v.s. (d,h) and (b) V5Q v.s. . The 
symbols denote the lattice QCD data, and the curves the theoretical form of the OGE plus multi-Y 
Ansatz. The lattice QCD results indicate the multi- Y-shaped flux-tube formation in the QQ-Q-QQ 
system. 



4. Heavy-Heavy-Light Quark Potential and Quark Motional Effect 

To see quark motional effects in baryons, we perform the first study of the hea yy- 
heavy-light quark (QQq) potential both in lattice QCE021 and in a quark modelM^ 

We calculate the QQq potential VQQq(-R) as the function of the distance R 
between the two heavy quarks in SU[3) quenched lattice QCD with (16 4 , (3 = 6.0) 
and 0(a)-improved Wilson fermionsffSl VQQ q (i?) is found to be well described with a 
Coulomb plus linear potential form, VQQ q (i?) = — ^rf" + u c sR + C c s, for R < 0.8fm. 
We find that, compared with the static three-quark case, the effective confining 
force (Jeff between the two heavy quarks is reduced in the QQq system by the qu ark 
motional effect, e.g., cr c ff — O.803Q for the constituent quark mass M q ~ 0.5GeVp^ 

We also calculate the QQq potential VQQ q (i?) in the quark potential model with 
the three-quark confinement potential obtained by lattice QCD.Q^The light-quark 
wave-function ip(r) distributes in the spatial region between the two heavy quarks 
QQ, as shown in Fig. 7. We find again reduction of the effective confining force <j c g 
between QQ in the QQq system, e.g., a Q g ~ O.80-3Q for M q ~ 0.5GeV.L^ 

Both in lattice QCD and in a quark model, the inter-two-quark confining force is 
effectively reduced by the motional effect of the remaining "3rd" quark in baryons. 
We conjecture that the effective reduction of the inter-two-quark confining force 
generally occurs also in light-quark baryons due to the 3rd quark motional effect. 
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Fig. 7. The light-quark spatial distribution \ip{r)\ 2 for M 9 =330MeV with _R=0.06fm (left), 
ij=1.2fm (center) and i?=2.4fm (right) in QQq systems. The brighter region has higher prob- 
ability, and the black circles denote the positions of the heavy quarks. 



5. Summary and Concluding Remarks 

We have performed the first detailed studies of various inter-quark potentials in 
SU(3) lattice QCD. The static three- and multi-quark potentials are well described 
by Y-Ansatz and the OGE Coulomb plus multi-Y type linear potential, respectively. 
This supports the flux-tube picture for baryons and multi-quark hadrons. We have 
also studied heavy-heavy-light (QQq) quark systems, and have found the quark 
motional effect to reduce effectively the inter-two-quark confining force in baryons. 
These lattice QCD results are useful to construct the QCD-bascd quark model. 
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